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The positron-neutrino correlation in the + — > + (3 decay 
of 32 Ar was measured at ISOLDE by analyzing the effect of 
lepton recoil on the shape of the narrow proton group follow- 
ing the superallowed decay. Our result is consistent with the 
Standard Model prediction. For vanishing Fierz interference 
we find a = 0.9989 ± 0.0052 ± 0.0036, which yields improved 
constraints on scalar weak interactions. 
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The e + -v correlation in + — > + (3 decay provides a 
robust signature of possible rare processes where parti- 
cles other than the usual W boson are exchanged. In 
the Standard Model and its left-right symmetric general- 
izations, + — > + decays produce the e + and v with op- 
posite chiralities. Angular momentum conservation then 
prevents relativistic leptons from being emitted in op- 
posite directions. On the other hand, a scalar interac- 
tion, which could arise from exchange of a leptoquark 
or a Higgs boson Q, will produce the e + and v with 
identical chiralities so that, in the relativistic limit, the 
leptons cannot be emitted in the same direction. In min- 
imal extensions of the Standard Model, Higgs couplings 
are too small to affect significantly the e-v correlation. 
However, supersymmetric theories with more than one 
charged Higgs doublet can accommodate sizable scalar 
couplings [jl] that are not ruled out by existing data Q . 

This paper reports new constraints on scalar weak in- 
teractions based on a precise measurement of the e-v cor- 
relation in the + — » + (3 + decay of 32 Ar, the only pure 
Fermi transition whose e-v correlation has been deter- 
mined with good precision. The simple spin structure 
of this decay permits tests for scalar interactions with- 
out complications from axial or tensor currents or from 
recoil-order effects. We consider a + — > + /3 + -decay 
Hamiltonian H 

(4>n1pp) (C S 1p v 1p e + C' s 1p v J 5 1p e ) , (1) 

which gives a decay rate 
^0--F(^^(l + a|co S ^ + 6-) 
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where E, p and m are the total energy, momentum and 
mass of the positron, E v the energy of the neutrino, Mj 
and Mf are the masses of the parent atom and daugh- 
ter nucleus, and F(Z,p) the Fermi function. We assume 
that the Standard Model provides an exact description 
of the IL ± exchange process and that Cy — C v || , but 
make no assumptions on the parity or time-reversal prop- 
erties of the scalar interaction. Then the e-v correlation 
coefficient, 
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are functions of C's and C" Q where 
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The C's will be complex if the scalar sector violates time- 
reversal invariance. 

The e-v correlation must be inferred from the recoil 
momentum of the daughter nucleus. This traditionally 
was measured directly, restricting high-precision results 
to a few favorable cases such as 6 He pj and neutron decay 
However, if the daughter nucleus is particle unsta- 
ble, the daughter momentum can be determined from its 
decay products. This allows one to study energetic light 
particles rather than slow heavy ions whose atomic and 
even chemical effects must be considered. Furthermore 
the transformation from the rest frame of the daughter to 
the lab amplifies the lepton recoil effects by a factor 2V/ v 
where V is the center-of-mass velocity of the light particle 
and v is the daughter's veleocity due to lepton recoil. Fi- 
nally, the time scale for particle decay is so short that the 
delayed particle is emitted before the recoiling daughter 
can slow down appreciably. Clifford et al. compared 
the energies of delayed en's detected in coincidence with 
P's emitted toward and away from the a, while Schardt 
and Riisager || studied the broadening of narrow delayed 
proton groups due to lepton recoil. The coincidence tech- 
nique allows one to measure a shift rather than a spread 
which is favorable on statistical grounds and is less sen- 
sitive to the response of the charged-particle detector. 
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On the other hand, the extracted value of the correla- 
tion coefficient is very sensitive to the energy and angle 
of the detected (3 particle. We adopted the singles tech- 
nique because it seemed difficult to determine the beta's 
kinematics with sufficient precision. 
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FIG. 1. Intrinsic shapes of the + — > + delayed proton 
group for a — +1, b = (heavy curve) and a — —1, b = 
(light curve). The daughter's 20 eV natural width is not vis- 
ible on this scale. 

The kinematics of 32 Ar superallowed decay deter- 
mine the daughter's velocity distribution and thereby the 
broadening of the delayed proton peak. The maximum 
kinetic energy of the recoiling 32 CI nucleus is 




2M, 



(6) 



where A is the difference in atomic masses of the par- 
ent and daughter states. The accepted value of A is 
poorly known because of the ±50 keV uncertainty in the 
32 Ar mass H. We use the isospin-multiplet mass equa- 
tion Jh}, M(T 3 ) = c + ciT 3 + c 2 T 3 2 , to obtain an im- 
proved value for A from the measured masses of all 5 
members of the A = 32, T = 2 multiplet. These masses, 
shown in Table El were obtained from the known ground- 
state masses and excitation energies jl^] of the isobars, 
except for 32 CI which we computed from our measured 
energy of the superallowed delayed-proton peak in 32 Ar 
decay, E LAB = 3349.9 ± 1.2 keV @ and the known pro- 
ton and 31 S masses. The isospin-multiplet mass equa- 
tion provides an excellent fit to the data and predicts 
that A = 3c 2 - ci = 6087.3 ± 2.2 keV @, implying 
T max w 638 eV and a maximum daughter velocity of 
v = 2.07 x 10~ 4 c. The 32 C1 daughter state has a width 
r w 20 eV (see below) so that in one mean life the 
daughter travels at most 2.1 x 10 _2 A before emitting 
the proton. The recoiling 32 CI therefore emitted the pro- 
ton while it was still traveling with the full velocity it 
received from lepton recoil. The intrinsic shape of the 
delayed proton peak (the shape for a counter with per- 
fect energy resolution) is shown in Fig. [l] for the limiting 
cases a = +1. 6 = and a = — 1, 6 = 0. 



We performed our experiment at ISOLDE. Beams of 
60 keV 32 Ar and 33 Ar ions from the General Purpose 
Separator were focused through a 4 mm diameter colli- 
mator and implanted in a 22.7 /ig/cm 2 carbon foil in- 
clined at 45° to the beam axis. Protons were detected in 
a pair of 9 mm x 9 mm PIN diode detectors collimated 
by 7.72 mm x 7.72 mm apertures located 1.6 cm from 
the beam axis. We eliminated possible uncertainties from 
beta summing effects by placing the detection apparatus 
inside a 3.5 T superconducting solenoid. The magnetic 
field prevented the betas from reaching the proton detec- 
tors (the highest energy betas from the + — * + decay 
had R c = 0.53 cm), but had little effect on the protons 
(the superallowed proton group had R c = 7.56 cm). 

The PIN diodes were maintained at — 1 1 C by thermo- 
electric elements that held the diode temperatures con- 
stant to ±0.02 C. The signals were amplified by preampli- 
fiers located immediately outside the vacuum chamber. 
The preamplifier housings were held at ±20 C by ther- 
moelectric devices that held the housing temperatures 
constant to ±0.01 C. Condensation of vacuum system 
contaminants on the detectors and stopper foil was min- 
imized by surrounding them with a copper shield cooled 
by a steady flow of liquid nitrogen. As an added pre- 
caution, the detectors were warmed to +27° C once each 
day to drive off any condensed material. The preamplifier 
signals were amplified and digitized by modules mounted 
in temperature-controlled crates and recorded in event- 
mode by a mini-computer. For each event we recorded 
the detector energy signals, the absolute time, the delay 
time after the arrival of a proton pulse, and the tem- 
peratures of the detectors, preamps, NIM crate, liquid 
nitrogen shroud, and the room. Our system gave excel- 
lent resolution; the pulser peaks for the two detectors had 
full-widths at half-maximum of 2.98 and 3.27 keV. 

Data were taken over a period of 12 days under several 
different conditions: with the stopper foil at 45°, 135°, 
225° and 315° with respect to the beam axis, and for two 
different beam tunes. These produced 6 different spectra 
for each of the 2 counters. We continually alternated 
between w 2 h long 32 Ar runs and 5-15 min long 33 Ar 
runs that provided energy calibrations for the 32 Ar data. 
The 32 Ar and 33 Ar beam intensities on target, 

TABLE I. Comparison of the measured mass excesses of 
the lowest T = 2 quintet in A = 32 to predictions of the 



Isobaric Multiplet Mass Equation [P(x 2 ,^) 


= 0.71]. 


isobar 
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M cxp (keV) a 


Mimme (keV) 


32 Si 


±2 


-24080.9 ± 2.2 


-24081.9 ± 1.4 


32 p 


+1 


-19232.88 ± 0.20 b 


-19232.9 ± 0.2 


32 g 





-13970.98 ±0.41 c 


-13971.1 ±0.4 


32 a 


-1 


-8296.9 ± 1.2 d 


-8296.6 ± 1.1 


32 Ar 


-2 


-2180 ± 50 


-2209.3 ±3.2 



a unless noted otherwise, ground state masses are from Ref. 

b E x = 5072.44 ± 0.06 keV from Ref. Ill 

°E X = 12045.0 ± 0.4 keV from Refs. @0 

d from delayed proton energy Jll| and masses of Ref. ||. 
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averaged over the entire run, were 94 and 3900 ions/s. 
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FIG. 2. Fit (upper panel) and residuals (lower panel) of the 
+ — ► + delayed proton peak. This spectrum, (the sum of 
detector 2 data in reflection geometry) contains roughly 1/4 of 
our data. The energy scale is 0.500 keV/channel. The pulser 
peak shows the electronic resolution. The Breit-Wigner tail 
from the 20 eV daughter width is visible on the high-energy 
side of the peak. 

We computed the intrinsic proton shapes using Monte 
Carlo routines to generate /3-decay events distributed ac- 
cording to Eq. U with 

P _ (-5m — E) — A — T (7\ 

1 + (pcos9 eL , - E)/Mi 

The Fermi function for a screened, finite-sized nu- 
clear charge was interpolated from Tables II and III of 
Ref. fTojl , and Gliick's Q order-a radiative correction 
to the energy distribution of recoiling 32 CI nuclei was ap- 
plied. A predicted (lj] 6.7 x 10~ 4 electron-capture branch 
was also included. Protons were ejected isotropically in 
the CI frame and deflected by the magnetic field; the 
mean energy losses of individual protons in the stopper 
foil and detector dead layer (roughly 1.5 and 1.8 keV, 
respectively) were computed. The stopper foil thickness 
was deduced from the energy loss of 3183 keV 148 Gd a's 
in the foil, while the implantation profile of Ar ions in the 
foil was computed with TRIM [Q. The detector dead 
layer thicknesses, 23.4 ±0.4 and 21.6 ±0.7 /xg/cm 2 , were 
measured by inserting into the apparatus a jig that al- 
lowed the 148 Gd source to be moved along an arc centered 
on one of the PIN detectors. The dead-layer loss varied 
as the secant of the angle while the energy loss in the 
source was essentially constant because the a source was 
always perpendicular to the line of sight to the detector. 

We evaluated the intrinsic proton shape at 961 points 
on a Cs-C'g grid. The intrinsic shapes plus a small, flat 
background were convoluted with a proton detector re- 
sponse function consisting of 2 low-energy exponential 
tails folded with a Gaussian. This functional form gave a 



good parameterization of a "first principles" calculation 
of the response function that included the effects of pulser 
resolution, the Fano factor, electronic and nuclear strag- 
gling in the stopper foil and detector dead layer, escape 
of Si X-rays, and energy loss to phonon excitations of the 
detector. We fitted our 6 pairs of delayed-proton spectra 
by varying the response function parameters (Gaussian 
width, tail lengths and fractional areas) to minimize x 2 
for each Cs, C' s point. Figure || displays the quality 
of the fits. Figure || shows that the extracted lineshape 
agreed well with the "first-principles" calculation. Fits 
were made for a series of values of Y, the natural width 
of the daughter. There was essentially no correlation be- 
tween r and Cs and C' s , and we found Y = 20 ± 10 eV. 
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FIG. 3. Comparison of the detector response function ex- 
tracted from the data in Fig. ^| to the "first-principles" calcu- 
lation described in the text. 



Figure [| shows constraints on Cs and C' s from this 
work and from Refs. |l9|~[23]|. The annular shapes of 
our constraints arise from the Fierz interference term in 
Eq. |^. Our Cs-C' s constraints may be parameterized as 

a = a/ (I + 0.19136) 
= 0.9989 ± 0.0052(stat.) ± 0.0036(syst.) 68% c.l. (8) 

where a and 6 are given in Eqs. ||, || with (m/p) = 0.21. 
Note that a, unlike a, does not have an upper bound of 
+1, so the range spanned by our experimental 2a error 
band lies entirely within the physical region. 

The systematic error included in Eq. || and Fig. || was 
evaluated by combining in quadrature the following ef- 
fects. We found the dependence of a on the exact val- 
ues of A and Q v , da/dA = -1.2 x 10~ 3 keV -1 and 
da/dQp = —0.9 x 10~ 3 keV -1 , by repeating the en- 
tire analysis with A and Q p changed by ±10 keV. The 
uncertainties, SA — ±2.2 keV and 5Q P — ±1.2 keV, 
gave a kinematic systematic error 5a = ±0.0032. p3 We 
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checked the dependence of a on the fitting regions of the 
proton spectra; a 28% variation in the width of the re- 
gion changed a by less than ±0.00055. We examined the 
dependence of our results on the form of the detector 
response function by reanalysing the data with a single- 
tail response function; by reanalyzing the data assuming 
that a weak Gamow- Teller peak lay under the tail of the 
32 Ar superallowed peak; and by simultaneously fitting 
the 33 Ar and 32 Ar superallowed peaks using a common 
response function. From these tests we inferred a line- 
shape systematic error of 5a = ±0.0016. 
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FIG. 4. 95% conf. limits on Cs and C's- Upper panel: 
time-reversal-even couplings. The annulus is from this work. 
The narrow diagonal band is from b(0 + — > + ) The 
broad diagonal band shows constraints from A, B, a, and 
ti/2 in n decay [^; the sausage-shaped area includes, in ad- 
dition, constraints from G( 14 0) and G( 10 C) |l|, fe( 22 Na) || 
and a( D He) g. Lower panel: time-reversal-odd couplings. 
The circles are from this work and correspond to Cs and C' s 
phases of ±90° , +45° and —45° . The shaded oval shows the 
constraint with no assumptions about this phase. The diago- 
nal band is from #( 19 Ne) @. 



coupling strength is M s = \C S \~ 1/2 M W > A.1M W . 
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